Resonating Valence Bond Mechanism of 
Impurity Band Superconductivity in Diamond 
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Superconductivity in an uncompensated boron doped diamond, a very recent observation, is strik- 
ingly close to an earlier observation of Anderson-Mott insulator to metal transition, prompting us 
to suggest an electron correlation driven superconductivity in an impurity band. Random coulomb 
potential remove a three fold orbital degeneracy of boron acceptor states, resulting in an effective 
single, narrow, tight binding and half filled band of holes. Singlet coupling between spins of neigh- 
boring neutral acceptors B° — B° is the seed of pairing. Across the insulator to metal transition, 
a small and equal fraction of charged B + and B~ states (free carriers) get spontaneously gener- 
ated and delocalize. Thereupon neutral singlets resonate and get charged resulting in a resonating 
valence bond (RVB) superconducting state. 

PACS numbers: 



Diamond has been reported^] to become a supercon- 
ductor upon high boron doping with a T c ss 4 K. This 
remarkable discovery of superconductivity in an uncom- 
pensated p-type semiconductor has possible implications 
for basic science and technology. As an example, shallow 
dopant states in Si and Ge play a key role in modern 
solid state electronics; from physics point of view they 
offer a testing ground for various ideas such as Anderson 
and Mott localization and their interplay. Boron doped 
diamond 0, being a relatively deep level acceptor with a 
hole binding energy of » 0.37 eV provides a new energy 
scale, new possibilities Q and a rich physics. Further in- 
vestigation of this system should be rewarding. 

The present problem of electronic phases of boron 
doped diamond, containing competing coulomb interac- 
tion, randomness and coupling to phonon, is like Si:P, 
quite a complex and challenging many body problem^. 
However, using a body of insights one has gained in the 
last couple of decades from theory and experiments one 
can suggest physically plausible mechanism and insights, 
which could be guiding hypothesis to understand the ob- 
served superconductivity in boron doped diamond (dia- 
mond:B). It is in this spirit, the present letter suggests 
a minimal model and mechanism for superconductivity. 
The dopant density seems 7] to be just above the critical 
density of the Anderson-Mott transition making electron 
correlation important. We suggest a resonating valence 
bond (RVB) mechanism of superconductivity in 
boron impurity band system. Diamond, a broad band 
insulator, where electron correlations are not important, 
should be viewed as offering an appropriate vacuum to 
boron (atom with an odd electron number) subsystem, 
where electron correlation driven Anderson-Mott insula- 
tor to superconductor transition and other rich correla- 
tion physics could take place(figure 1). 
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FIG. 1: Hole density of states (schematic) in boron doped 
diamond, an uncompensated p-type semiconductor. Holes of 
acceptors form a strongly correlated and effective single im- 
purity band at half filling. Anderson-Mott insulator to super- 
conductor transition is suggested to take place in the impurity 
band as we increase boron density (figure 2). 



The observed superconductivity in diamond:B is of 
type-II with a T c w 4 K, H C 2 > 3.5 T and coherence 
length £gl ^100 Au. The dopant boron concentration 
is around 4.5 x 10 21 /cm 3 . 

Transport properties of diamond:B has been studied 
by various authors and the issue of metal insulator tran- 
sition in the impurity band discussed 0, E|- The only 
report to my knowledge, which makes a low tempera- 
ture measurement of the critical boron concentration by 
a scaling analysis is by Tshepe et al.[7j. They find a 
critical doping concentration n c s» 4 x 10 21 /cm 3 , strik- 
ingly close to the boron concentration ~ 4.5 x 10 21 /cm 3 
used in the new superconductor It is this closeness 
which prompted us to examine the role of electron corre- 
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lation and suggest a Anderson-Mott insulator to RVB su- 
perconductivity transition, similar to the one advocated 
recently in the context of (non-random) Mott insu- 
lator to superconductor transition. Importance of strong 
correlation is further enforced by a nearly temperature 
independent and a very large value of normal state resis- 
tivity of the (polycrystalline) superconducting sample^], 
p ~ 9 mfl cm, which gives a mean free path £q ~ 0.1 Au, 
much less than the average separation between neighbor- 
ing acceptor sites. Resistivity in single crystal or low 
frequency conductivity measurements in polycrystalline 
sample will be welcome to substantiate this point. 

Diamond, an excellent insulator has a very broad 
band and a wide gap ~ 5.6 eV. The top of the va- 
lence band is three fold degenerate with effective masses, 
2.12 m e , 0.7 m e and 1.06 m e . Boron, which has one 
less electron compared to a carbon, becomes an accep- 
tor, when it substitutes a carbon atom in diamond[2j. 
The acceptor states are threefold degenerate with a hole 
binding energy of Eb ~ 0.37 eV. Spin-orbit coupling 
partly removes this degeneracy by a marginal 6 meV. 

The relative dielectric constant of diamond is eo ~ 5.6. 
Effective mass theory estimate of acceptor states gives a 
hydrogenic impurity state with an effective Bohr radius 
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of a* B = 2 eoE B ~ ^ Au, of the 'envelope function'. When 
this is used in conjunction with Mott's criterion of Mott 
insulator metal transition, n c a 3 rj 0.25, one obtains (ll) 
a critical doping, n c « 8.0 x 10 21 /cm 3 . Experimentally 
observed critical concentration of reference [7] is of the 
same order as this rough estimate. 

Keeping the above in mind, we build a simple tight 
binding model for the impurity band. The three fold 
degenerate acceptor wave functions (transforming as 
t2g) have different spatial density distribution \ipa(Y — 
Ri)| 2 ,for a = 1,2 and 3. Here Rj is the site of a boron 
atom. Hence the total electrostatic potential a hole in an 
acceptor state feels from neighboring dopants depends 
on the orbital it occupies. In other words, the three fold 
degeneracy of an acceptor state is in general lifted by 
the randomness. Simple estimate shows that the amount 
of lifting, for the boron doping of interest, is large com- 
pared to the impurity band width ~ 0.2 eV. For low 
energy physics of interested to us, only the lowest energy 
acceptor states are important. This is also consistent 
with general experience in Si:P, where the degeneracy of 
donor states has no substantial role and a non-degenerate 
donor state theory seem to work wellfl^. 

Once we identify a lowest energy acceptor state for 
each boron, a one band tight binding model Hamiltonian 
follows: 

i 



Here c\ a is the hole creation operators at the lowest en- 
ergy acceptor state of boron at site i with an energy ej 
and tij are the hopping matrix elements; rij = + n^. 
The parameter f7, is the Hubbard repulsion term of an 
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acceptor state centered at site i and Wa ~ — %— is the 
diagonal coulomb matrix clement between two acceptor 
states separated by a distance Rij . The total number of 
holes is the same as the number of boron atoms; we have 
a half filled band of interacting holes. 

On the insulating side of the above Hamiltonian, the 
low energy sector is the spin sector, which can be seen 
clearly in the limit U » t as a random antiferromagnetic 
Heisenberg Hamiltonian with superexchange interaction 
Jij: 

ij 

This is well established in the context of Si:P for example 
leading to notions of hierarchy of singlet spin coupling, 
valence bond localization, valence bond glass etc., both 
experimentally and theoretically^^. We do not think 
diamond:B is fundamentally different from Si:P in the 
Mott insulating side. 

The issue for us is the conducting side close to Mott 
transition point. In the present paper, based on our re- 
cent study, we offer a new insight. Conventionally the 
conducting side is thought of as a disordered fermi liquid 
and with regions of certain local moment character [14j. 
We have recently suggested |9( that a corresponding con- 
ducting state close to the Mott transition point in the 
non-random case is well thought of as a 'self doped Mott 
insulator', a projected metal. In a self doped Mott insu- 
lator a small and equal density of free positive (B + ) 
and negative (B~) carriers are self consistently gener- 
ated across the Mott transition, out of the neutral states 
B° , in the process of minimizing the free energy, partic- 
ularly the long range coulomb interaction part. Further, 
antiferromagnetic superexchange, that is characteristic 
of a Mott insulator survives in the conducting state as 
well, the same way superexchange survives in the C1LO2 
planes of cuprates in the presence of doped holes. That 
is, in the above conducting state, in addition to 'virtual' 
double occupancy and empty sites, which are responsible 
for generation of superexchange, a small and equal den- 
sity of 'real' and delocalized double occupancy and empty 
sites are maintained. Our workQ unified RVB mecha- 
nism of superconductivity in hole doped cuprates with 
that in organic superconductors such as 2D ET salts and 
also predicted new systems. 

An effective Hamiltonian of the above conducting 
state is the 2 species random t-J Hamiltonian that is 
obtained^ by a superexchange perturbation theory for 
the present situation of half filled band of disordered in- 
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teracting electron Hamiltonian (equation 1) 

H U = ~ Y % Prf C ta C ja P d - Y % Pe C L C j« P e + h - c - 
ij ij 

- Y J v( Si ■ s j ~ \ n i n j) + Y eic ° c ° ' ( 3 ) 

ij i 

operating in a subspace that contains a fixed number xN 
of doubly occupied B~ and xN empty B + states. Here TV 
is the total number of electrons, which is the same as the 
number of lattice sites. The projection operators Pa and 
P e allow for the hopping of charged B~ and B + states 
respectively, maintaining (1 — 2x)N singly occupied neu- 
tral states B° in a dynamical fashion. Notice that the 
long range coulomb interaction has disappeared in the 
above Hamiltonian because of screening in the metallic 
state; More importantly it has left its effect parametri- 
cally through the generation of a finite density 2x of free 
carriers by the process of self doping. 

Our 2 species random t-J model adapted to the self 
doped Mott insulator has a more transparent form in the 
slave boson representation cJ CT = s\ a d i -\- as is e\. Here the 
charge-ons cij's and ej's are hard core bosons that create 
doubly occupied sites (B~ of charge e~) and empty sites 
(B + of charge e + ) respectively. The fermionic spinon 
operators sj CT 's create singly occupied sites (B° , charge 
neutral) with a spin projection a. The local constraint, 

<H^i + e l e i + J2a s \a s ia — ^ keeps us in the right Hilbert 
space. 

In the slave boson representation our t-J model takes 
a suggestive form: 

H t j = - Y *« ( e K- Y s tv - d Uj Y s L s j*) + h - c - 

ij a a 

- E^lA.+E^H-ek) (4) 

ij i 

where b\. = ^(4t s ]j. _ -sj^ s 1t ^ is a spin sm §l et spinon 
pair creation operator at the bond ij. 

At a formal level, a slave boson analysis is readily per- 
formed for our random interacting electron model. We 
suggest two different phases, a valence bond glass (VBG) 
and RVB superconducting phase, shown in figure 2. 

In the Anderson-Mott insulator region we have a VBG 
phase. This is brought out well[13| both in experiments 
and theory in the case of Si:P. This state contains only 
real 'neutral' 73° states and no real B + and B~ charged 
states. A good variational wave function to describe such 
a valence bond glass state is 



\VBG)=P(Y<t>VBG^j)bl 1 ) N \0), 



ij 



(5) 



that describes condensation of neutral valence bonds in a 
specific valence bond pattern given by the pair function 
0vbg(*j j)- While one can formally write down a self 
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FIG. 2: Schematic Phase Diagram as a function of dopant 
density in Diamond:B, an uncompensated case. 



consistent gap equation for the pair function </>vbg(*jJ) 
of our random correlated spin-half system, one does not 
attempt to solve it but assumes that a solution exists. 

In Si:P, the hierarchically organized valence bond 
couplings[13| in a VBG leads to a 'pseudo gap' that gives 
a magnetic susceptibility that vanishes as Xspin ~ T \- a 
with a > 0. The spin contribution to specific heat also 
gets a similar power law correction. 

As mentioned earlier, across the Mott transition a 
small and equal density of positive and negative carri- 
ers are spontaneously generated out of the insulating va- 
lence bond glass state. The derealization of charged 
carriers causes charging and resonance of the frozen va- 
lence bonds resulting in an RVB superconducting state. 
This state is best described by a generalized RVB vari- 
ational wave function, inferred from a slave boson mean 
field analysis: 



\RVB) 



p (eir N (4r N (bi)^ N \o) (6) 



Here e£ = 7^ d>*(i)e] and d\ = Y^iKi^A represent 
'bose condensation' 1 1 5| of the B + and B~ in two different 
'extended states' </> M («) and <j> v (i). As the hopping ma- 
trix elements of B + and B~ have opposite sign (equation 
4), they condense in general in different extended states. 
The operator bo = yj 4 - <psc(i,j)b\j represents the con- 
densation of valence bond pairs B° — B° in an extended 
state represented by the pair function (f>sc{i,j)- And P 
is a projection operator that ensures the presence of only 
one of B° or B + or B~ acceptor states in any boron site 
in the many body wave function, by obeying the local 
constraint, d\d. t + e\e t + J2a 4* s i* = l - 

To understand the above variational wave function 
(equation 6), we wish to state that in a non random case 
of simple cubic lattice for example, the holons and dou- 
blons respectively condense at wave vectors (0,0,0) and 
(7r,7r,7r), in view of the different signs of the holon and 
doublon hopping matrix elements in equation 4. Also 
note that in the RVB theory, the holon and doublon con- 
densation, an apparent charge 'e' condensation [l5j| is ac- 



4 



tually a book keeping device 16] for the charge '2e' con- 
densation of physical electron pairs (valence bond pairs) . 

As mentioned earlier, in view of the random charac- 
ter of our Hamiltonian, one can only make some exis- 
tence type of statement of the functions <j>fi(i), 4> v (i) and 
4>sc(hj) from plausibility arguments. For the same rea- 
son it is difficult to get quantitative estimate of the su- 
perconducting T c , from a formally 'exact' gap equation, 
even within the slave boson mean field analysis. How- 
ever, we can get a very rough order of magnitude estimate 
of the superconducting T c using certain 'typical' values 
of the impurity band parameters using a standard RVB 
mean field expression, ksT c « ^-e~~ . If we assume 
an impurity band width W around 0.2 eV and superex- 
change J around 0.04 eV, we get a T c in the range 1 to 
10 K. It is important to remark that what is seen ex- 
perimentally in diamond:B is likely to be some kind of 
lower bound for a larger intrinsic superconducting T c , 
as inhomogeneities in the real system are likely to affect 
T c considerably. 

What is the order parameter symmetry ? Super- 
conductivity is in a disordered lattice of boron accep- 
tors. Unlike cuprates, where a d x 2_ y 2-w&ve symmetry is 
compatible with an underlying square lattice symmetry, 
higher angular momentum states have no compatibility 
with the underlying random lattice. Hence, diamond:P 
is likely to go into an extended s-wave superconduct- 
ing state, there by also respecting Anderson's theorem 
on dirty superconductors. 

The conducting side of the Anderson-Mott transition, 
in view of residual unscreened short range coulomb re- 
pulsion could create a charge density glass of the spon- 
taneously created B + and B~ carriers as well. This will 
be a competitor to superconductivity. 

If our correlation mechanism is at work, the normal 
state should be an anomalous one. Disorder should lead 
to certain soft and localized spinon and chargeon quasi 
particle states. This is likely to lead to normal state 
anomalies different from cuprates. 

Our current proposals, including mechanism of super- 
conductivity, as it stands is very heuristic and based on 
limited available experimental results. To make further 
progress and check the validity and consequences of the 
present proposal, it is important to perform certain ex- 
periments: i) precise measurement of the critical con- 
centration n c of the Anderson-Mott insulator to metal 
transition ii) map out the superconductivity phase di- 
agram by changing the boron concentration iii) look for 
anomalous T a type power law signals in specific heat and 
magnetic susceptibility in the insulating phase, charac- 
teristic of the valence bond glass phase and hierarchical 
singlet coupling, iv) 'spin gap' type of behaviour in the 
normal state on the conducting side, looking for onset 
of spinon pairing v) normal state transport anomalies on 
the conducting side and vi) look for superconductivity 



in the Mott insulator phase by small amount of compen- 
sation, say by nitrogen type donor impurities; this will 
be an external doping of the boron impurity band rather 
than self-doping. 

Our present proposals raise some questions and sug- 
gests applicability to other systems. Why superconduc- 
tivity is absent in Si:P and related systems ? Apart from 
the fact that the energy scales are low, the impurity state 
orbital degeneracy is high and spin-orbit coupling is rel- 
atively high. Hund coupling and reduced band filling re- 
sulting from multiple bands is likely to diminish possibil- 
ity of singlet superconductivity. It will be also interesting 
to look for new systems and also find out if electron corre- 
lation play any important role in estabilishing supercon- 
ductivity in known impurity band superconductors [1J] 
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